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Based on the developed cross-correlation theory for remote location of neutrino source with two- 
detector setup for neutrino detection the modification of arrangement of the OPERA experiment 
is suggested. Within the framework of computing experiment based on the OPERA experimental 
data we show that the use of this theory makes it possible not only to determine with high accuracy 
the delay time between neutrino signals but to eliminate the errors of blind analysis by which all 
necessary time corrections for determination of signal "technologically unremovable" delay time 
between CERN and GSL are performed. 



INTRODUCTION 

Based on data processing over the period 2009-2011 the 
OPERA collaboration has recently announced the results 
about possible evidence for superluminal propagation of 
neutrinos [1]. 

As the recent developments show [2], the most weak 
point of the OPERA data interpretation [1] associated 
with the superluminal neutrino detection, is so-called 
blind analysis by which all necessary time corrections 
are performed and thereby the total signal "technolog- 
ically irremovable" delay time between CERN and GSL 
is determined. 

In this paper we propose the cross-correlation theory 
for remote location of neutrino source (e.g. CERN Target 
chamber) based on the two-detector (e.g. GSL and New 
Detector GSL-clone) setup for neutrino detection, which 
makes it possible not only to determine the delay time 
between neutrino signals, but to eliminate the errors of 
blind analysis as well. 



CROSS-CORRELATION THEORY OF 
NEUTRINO SOURCE LOCATION 

The proposed method for remote location of neutrino 
source is based on an intensity interference correlation of 
neutrino signal sequences from the same source measured 
by two or more spaced-apart detectors. 

First of all, let us say a few words about the cause 
of hidden but really existent interference of neutrino in- 
tensities in the OPERA-experiment. It is known, that 
within the framework of this experiment muon neutri- 
nos are produced due to pion and kaon decays, which, 
in its turn, are produced in proton-proton collisions at 
a graphitic target in CERN. It is obvious that in the 
OPERA experiment interference correlations of identical 
particles [3], which a priori take place in hadron-hadron 
interactions [4] , are " genetically inherited" by muon neu- 
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FIG. 1. Schematics of neutrino source passive location. See 
explanation in the text. 



trinos. Just these inherited properties of neutrino became 
the basis for the proposed theory. 

Let us consider the idea of interference location [5, 6] 
based on the simple scheme of spaced-apart neutrino de- 
tectors (Fig. 1). Note that by "neutrino detector" we 
mean a system, in which a number of measurable nuclear 
physical transformations, causally related to the detected 
neutrino, takes place. 

Within the framework of such statement of the prob- 
lem we are interested in an interference correlation of 
neutrino intensities generated by the wave functions ipi 
and ip2 , which describe the neutrino motion from a source 
to the detector 1 and detector 2, respectively (Fig. 1). In 
our case, a radially divergent wave from a source at the 
distance R, which significantly exceeds the source size, is 
described by the normalized spherical wave function. If 
the emission process lasts long enough, this function has 
the form of monochromatic spherical wave 
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where f{6) is the angular distribution amphtude of the 
radiation from source which depends on the angle 6 of 
divergent spherical wave; E and p is the neutrino energy 
and momentum; if is the initial phase of wave function; 
R is the distance between neutrino source and detector. 

Since the time of each neutrino production is sufh- 
ciently short, the wave function differs from a monochro- 
matic wave and looks like 
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Then the pulse sequence of neutrinos coming one by 
one to the detector 1 corresponds to the following total 
wave function 
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and the signal intensity measured by the detector 1 is 
determined by the correlation 



Since the sequence of randomly distributed pulses mea- 
sured by each of the detectors is the Gaussian process, 
Eq. (7) can be transformed according to a known rule 
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As a result, the cross-correlation function for signal 
intensities (7) measured by the two identical detectors 
has the following form 
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For the Gauss stationary processes the reciprocal 
moments {^i{t)£,i{t + t{, {Ci {t)£,* {t + t{ of the random 
events Ci{t) ^nd ^j^i{t + r) are zero. As a result, we 
have 
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where \^i{t)\ is the neutrino flux intensity in the detec- 
tor 1; El is the registration efficiency of neutrino events 
by the detector 1. 

By analogy, the neutrino pulse sequence in the detec- 
tor 2 corresponds to the following total wave function 



^2{t + T)^J2 ^2 {R2,t- t2n + t) 



E4?^^(^2,i-t2„+T)e^^^ (5) 



47ri?2 

where r is the additional delay of neutrino on the way 
from a source to the detector 2 with respect to the neu- 
trino traveling time from a source to the detector 1. The 
signal intensity measured by the detector 2 is determined 
by expression 
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At first let us consider the operation of correlation 
scheme disregarding the influence of filters 1 and 2 and 
delay line shown in Fig. 1. The cross-correlation function 
of signal intensities measured by detectors looks like 

Kj (r) = ( Ji (t) J2 {t + t)) - ( Ji {t)) (J2 [t + t)) = 
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From (10) it follows that the intensity cross-correlation 
function is equal to the intensity of the amplitude cross- 
correlation function ^(i). For the considered stationary 
processes, representing the random pulse sequences in the 
formulas (4) and (6), the central moment {^iit)^2it + ''')( 
is calculated in the standard way 
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where (N^) is the average number of neutrinos from a 
source during one measuring. 

Since the average number of neutrino events measured 
by any of the detectors per unit time is 
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the correlation function of neutrino signal intensity (9) 
will have the following dimensionless form: 
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Kj (r) = (m) (712) 
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Now let us analyze the operation features of considered 
scheme (Fig. 1) and possible ways of its optimization. We 
expand the random variable F{t) into the Fourier integral 
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ubstituting the expansion (14) into (13), we obtain 
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Then using the integral property 
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we transform the correlation function (15) into the form 



Kj (t) = (m) (n2) {2Tif X 

00 

j \F{Ri,u)F* {R2,uj)\\-'"^duj 



(17) 



This result is, in fact, a generalization of the Wiener- 
Khinchin theorem [7], which describes the time depen- 
dence of amplitude correlation function of random pro- 
cess, for the case of intensities interference (i.e., the 
quadratic function of amplitudes) of the same random 
process. 

Eq. (17) corresponds to the ideal correlator, which does 
not distort the processed signals F{Ri,uj) and F2 (i?2, w) 
from a detector. In the real experiment such distortions 
can be taken into account by introducing of the correlator 
transfer function G{uj) which yields the equation (17) in 
the form 
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Within the framework of OPERA experiment, the spec- 
trum width Suj w 27r/5t of the impulse signal F{R,t), 
formed at the nuclear particles detector output, may 
be considered as exceeding the electronic post-processing 
circuit transmission bandwidth Aw substantially. In this 
case one may assume that the correlator transfer func- 
tion corresponds to a filter with a rectangular frequency 
window Aw Slu and that the spectral density of each 
signal Fi 2{Ri^2, ^) is constant within this frequency win- 
dow and equals to -Fi,2(-Ri,2, 0). As a result, from (18) 
we obtain 
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To enhance the usability of this formula and precision 
of measurements we supplement the signal processing 
circuit with two specialized controlled frequency filters 
with the guaranteed rectangular pass band Awi < Aw 
(Fig. 1). Adding the delay line to one of the neutrino 
intensity interferometer arms (Fig. 1) will result in the 
expression for the interference correlation function (and 
the correlator output signal) in the form 
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The magnitude Kj (r — Ai) tops at the optimal pulse- 
delay time Ai = t. Thus, obtaining the maximum value 
of the correlator output signal Kj{T—llt), we experimen- 
tally determine the optimum pulse-delay time t. Below 
we consider this problem within the framework of com- 
puting experiment based on the OPERA-data. 



THE OPERA-DATA AND RESULTS OF 
COMPUTING EXPERIMENT 

At first, we consider the computational results of com- 
puting experiment on the statistical simulation of neu- 
trino measurement by two detectors located at different 
distances from a neutrino source. It is obvious that a 
statistics of neutrino detection must obey the cascade 
statistics of the Poisson compound distributions [8], so 
far as it is known, that the statistics of hadron multiple 
production in pp- and pp-collisions on the energy inter- 
val in c.m.s. -/s = 20 — 1800 GeV is well described by 
the known cascade distribution p(n) of the Neyman type, 
which belongs to the class of Poisson compound distri- 
butions [7]. We apply this distribution without the loss 
of generality in the simplest (for (n) 3> 1) recurrent form 
[7, 9-11]: 
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where (A) is the average number of proton-target cohi- 
sions; (e) is the average number of secondary particles 
(pious and kaons) per one primary coUision; p^ is the 
BernoulH probabihty of neutrino production due to sec- 
ondary particle decay. 

The average number of detected neutrino and variance 
have the form 
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The experimental distribution of neutrino events ob- 
tained by the OPERA-data for the first extraction [1] 
is shown in Fig. 2a. In spite of small length of neu- 
trino event sampling, this distribution is well described 
by the Neyman distribution with the same parameters, 
i.e., (n) = 113 and var{n) = 144. 

The computing experiment on the statistical simula- 
tion of neutrino measurement by two detectors located 
at different distances from a neutrino source (Fig. 1) was 
performed in the following way. Statistics of the neu- 
trino event number is described by the Neyman distri- 
bution (21) and is set by the random number generator 
with the average (n) = 113 and variance var{n) = 144 
(see Fig. 2). The average number of neutrinos from a 
source during one measurement is set so that the aver- 
age number of neutrino events (712) = 113 (Fig. 2) was 
registered by the detector 2 (Fig. 1). We considered two 
locations of the detector 1 (Fig. 1). In the first case the 
location of the detector 1 and neutrino source coincided 
(Ati2 = 100 a.u.), while in the second case the detector 1 
was located at the distance from a source, which corre- 
sponded to the propagation time Ati2 = 50 a.u. between 
the first and second detectors. The computational re- 
sults for the cross-correlation function of neutrino event 
intensity for two variants of the first detector location are 
presented in Fig. 2b. 

It is obvious also that in the general case of two inde- 
pendent detectors (for example, GSL and New Detector 
GSL-clone), which are located at the different distances 
from a neutrino source, the neutrino velocity is deter- 
mined as 
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where t,^ - is the delay time between the neutrino signals; 
Li and L2 are the distances between a neutrino source 
and detectors, respectively. 
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FIG. 2. (a) Comparison of the experimental distribution of 
neutrino events for the first extraction [1] with the average 
(n) = 113 and variance var{n) = 144 (blue histogram) with 
the Gauss distribution (green fill area) and Neyman distribu- 
tion (red histogram) with the same parameters; (b) the cross- 
correlation function (17) for the two samplings of neutrino 
events measured by two detectors spaced on the conditional 
time T — 50 a.u. and r — 100 a.u.. See explanation in the 
text. 



Let us recall that in our case the value are calculated 
by the expression for cross-correlation function (20). 

In conclusion, we give the quantitative evaluation of 
cross-correlation function for the OPERA-experiment. 
For this purpose we suppose that the proton extrac- 
tion waveform (Fig. 3a) measured by a fast Beam Cur- 
rent Transformer (BCT) detector [1] and transmitted by 
telecommunications to GSL in the form of the proton 
probability density function emulates some neutrino in- 
teraction time distribution "measured" in the place of 
the OPERA-detector location. By virtue of the linear de- 
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pendence between the proton beam intensity and average 
number of secondary hadrons (and, consequently, neutri- 
nos) this supposition is sufficiently admissible. With con- 
sideration of these assumptions we will consider bellow 
that the signal of "the proton probability density func- 
tion" travels between CERN and GSL with the speed of 
light, while the neutrino velocity is unknown. 

Thus, based on data shown in Fig. 3a we consider 
the problem of determining the delay time of one signal 
(the proton probability density function) with respect to 
other (the neutrino interaction time distribution) , which 
are measured in the place of the OPERA detector lo- 
cation. Computational results of the cross-correlation 
function of neutrino events intensity for the first extrac- 
tion are presented in Fig. 3b. One can see that with 
consideration of the BCT detector resolution {5t « 5 ns, 
Su! « 1.25-10^ Hz) our result (r = 1050 ns) is in complete 
agreement with the OPERA result {dti = 1043.4 ns). 



CONCLUSIONS 

Based on the developed cross-correlation theory for re- 
mote location of neutrino source we offer the new ar- 
rangement of experiment, by which the precision deter- 
mination of the neutrino velocity is possible. It is shown 
that the two-detector setup for neutrino detection makes 
it possible not only to determine with high accuracy the 
delay time between neutrino signals, but to eliminate the 
influence of blind analysis errors [1], accompanying the 
procedure of performing the necessary time corrections 
for determination of total signal "technologically unre- 
movable" delay time between CERN and GSL. 

It is necessary to note also a very substantial fact. Ap- 
plication of the cross-correlation method makes it pos- 
sible to measure with high accuracy the signal time de- 
lay without using the proton pulses with steep and sharp 
fronts. In other words, determination of cross-correlation 
function does not require the signal special preparation 
by the BCT-detector and allows one to use the arbitrary 
proton pulses with arbitrary "diffuse" fronts. This, in its 
turn, makes it possible to increase the statistics of neu- 
trino events and not to reject the results of trial starts 
even. 
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the Gran Sasso Underground Laboratory, where the idea 
of this paper arose. 
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